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Flavor SU(3) symmetry, including 30% first order SU(3) breaking, has been 
shown to describe adequately a vast amount of data for charmed meson de¬ 
cays to two pseudoscalar mesons and to a vector and a pseudoscalar meson. We 
review a recent dramatic progress achieved by applying a high order perturba¬ 
tion expansion in flavor SU(3) breaking and treating carefully isospin breaking. 
We identify a class of U-spin related D° decays to pairs involving charged pseu¬ 
doscalar or vector mesons, for which high-precision nonlinear amplitude relations 
are predicted. Symmetry breaking terms affecting these relations are fourth order 
U-spin breaking, and terms which are first order in isospin breaking and second 
order in U-spin breaking. The predicted relations are shown to hold within ex¬ 
perimental errors at a precision varying between 10~ 3 and 10“ 4 , in agreement 
with estimates of high order terms. Improved branching ratio measurements for 
D° —> K + p~, K* + n~ and for decay modes involving three other pairs of charged 
pseudoscalar and vector mesons could further sharpen two of these precision tests. 
We also obtain amplitude relations for D° decays to pairs of neutral pseudoscalar 
mesons, and relations for rate asymmetries between decays involving Kg and Kg, 
which hold up to second order U-spin breaking at a level of several percent. 


I Introduction 

A useful tool for studying hadronic decay amplitudes of charmed mesons is approximate 
flavor SU(3) symmetry. First order symmetry breaking corrections in amplitudes, due to 
the light quark mass term in the Standard Model Lagrangian, are expected to be of order 
(r», s - Tn u ,d) /Aqcd ~ fx/U - 1 ~ 0.2 - 0.3. 

In the seventies, shortly after the discovery of charm, SU(3) group theory has been used to 
obtain amplitude relations for charmed meson decays into pairs involving two pseudoscalar 
mesons or a pseudoscalar and a vector meson |1]. This study was extended in the early 
nineties to include numerous first order SU(3) breaking terms one of which was assumed 
to dominate over the others El El- A diagrammtic approach [3], equivalent to SU(3) group 
theory, has been developed and applied to fit data for branching ratios as they have been 
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accumulated [3]. Assumptions made in these studies about SU (3) breaking were often model- 
dependent. Other studies of these decay processes went much beyond SU(3) by assuming 
factorization of hadronic amplitudes [6]. A recent SU(3) fit to current data of charmed 
meson decays to two pseudoscalar mesons worked reasonably well when including two of the 
numerous first order SU(3) breaking terms (of order 30%) available in a group theoretical 
approach [7]. 

While flavor SU(3), including 30% first order symmetry breaking, has been shown to 
describe adequately a vast amount of data of hadronic decays of charmed mesons, this 
has not provided a precision test. For this matter one would hope to develop a perturbative 
expansion in reasonably small SU(3) breaking parameters, in which only high order symmetry 
breaking terms survive certain relations among amplitudes. A small step in this direction 
was made in Ref. [8], searching for linear relations among amplitudes for two-body and quasi 
two-body charmed meson decays in which first order SU(3) breaking terms cancel. Testing 
these relations, expected to hold within several percent, requires in most cases measuring 
relative strong phases between amplitudes which is a highly challenging task. 

The purpose of this paper is to review, expand and improve new results obtained in recent 
work published in two short reports 0 ED] applying a high order perturbation expansion 
in flavor SU(3) breaking and treating carefully isospin breaking. We will identify a class 
of D° decays, for which high-precision nonlinear relations among magnitudes of amplitudes 
hold. The lowest order symmetry breaking terms affecting these relations will be shown to 
be fourth order SU(3) breaking terms, and terms which are first order in isospin breaking 
and second order in SU(3) breaking [IT]. 

One major motivation for this work is searching for signals of new physics. Very precise 
relations as discussed here, which would fail at some high order flavor symmetry breaking, 
could provide such signatures. For great convenience we will use U-spiu, an SU(2) subgroup 
of flavor SU(3), rather than applying the full SU(3) group structure. Our U-spiu expansion 
is, of course, consistent with SU(3) expansions in Refs. 0 [7, S]. 

In Section [TT] we identify sets of four two-particle final states in D° decays, each consisting 
of pairs involving charged pseudoscalar and vector mesons. In a given set two of these states 
and a linear combination of the other two form a U-spin triplet, playing an important role 
in D°-D° mixing. Our discussion in Sections [111] - IVIII studies in detail one of these sets 
denoted by D° —> P + P~ consisting of D° —> n + K ~, K + n~, K + K~, 7r + 7r - . We derive U-spin 
symmetry relations for these processes in Section IIII1 and study first order and arbitrary 
order U-spin breaking corrections in Sections IIVI and [Vj respectively. In Section [VI] we 
obtain a high-precision relation, obeyed up to fourth order U-spin breaking, for ratios of 
amplitudes of the above four processes. Section IVIII investigates first order isospin breaking 
terms occurring in this relation, showing that they are also suppressed by factors associated 
with second order U-spin breaking. Sections IVIII! IIXI and |X] discuss experimental tests in 
D° —> P + P ~, D° —> V + P~ and D° —y P + V~, respectively, where V ± denote charged vector 
mesons. In Section|Xl]we study D° decays into pairs of neutral pseudoscalar mesons, deriving 
amplitude relations and Kg — K { } rate asymmetry relations which hold up to second order 
U-spin breaking. A short conclusion is given in Section Kill 
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II D° decays to two-body U — 1 states 


An SU(2) subgroup of flavor SU(3) that is useful for studying charmed mesons is U-spin [12]. 
The quark pair (d, s) behaves like a doublet under this group while the u quark and the 
heavier c, b and t quarks are singlets. U-spin symmetry leads to interesting relations among 
amplitudes of hadronic D decays nang. It also implies the vanishing of D°-D° mixing up 
to second order U-spin breaking [15], which underlies the vanishing of D°-D° mixing within 
full flavor SU(3) |16j . This behavior of D°-D° mixing under U-spin has been shown to 
follow from a cancellation up to second order U-spin breaking of mixing contributions from 
intermediate U-spin triplet states. The high order U-spin breaking perturbation expansion 
that will be studied in this paper works well, as we will show, for these two-body or quasi 
two-body D° decays to U — 1 states. 

One class of U — 1 two-body states involves pairs of opposite charge pseudoscalar or 
vector mesons. Our study will focus on these final states for which a high-order U-spin 
breaking expansion is applicable. Another class of processes involves decay into pairs of 
neutral mesons. In this case two-body final states do not have well-defined values of U-spin. 
Instead, linear superpositions of final states have U — 1. Consequently in these decays a 
U-spin expansion works well for certain linear combinations of decay amplitudes. 

We start by classifying single meson states of positive or negative charge as doublets of 
U-spin. Since a pair of d and s quark and their antiquarks form two U-spin doublets, ( d , s ) 
and (s, —d), one has two doublets of pseudoscalar mesons 



and two doublets of vector mesons 



( 1 ) 

( 2 ) 


One can then form two-particle states of charge zero in four different forms, P + P , V + P , 
P + V~ and V+V~. 

In the next several sections we will study the four processes D° —> P + P~, P + P~ = 
ti + K~, K + 7i~, K + K~, for which most precise data exist. This discussion is also 

applicable to the other three sets of processes, D° — > V + P~, P + V~ and D° — > V + V~. 
For D° -> V + V~ one may, in principle, treat separately S, P and D-wave amplitudes, or 
amplitudes for longitudinal polarizations, and for tranverse polarizations which are mutually 
parallel and perpendicular to each other. We will not discuss further these latter challenging 
decay modes for which no branching ratios have been measured di- 


iii U-spin symmetry limit 

As a starting point we derive amplitude relations in the U-spin symmetry limit. The four 
possible two-particle states |P + P~) can be written in the form of three U-spin triplet states 
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and one singlet state: 


\K + n~) = |1,1), |t t+K~) = -|1, -1), ^\ K+K ~ ~ ^ + 0 = |1, 0), (3) 

-±=\K + K- + n + n-) = \0,0) . (4) 

The charm-changing weak Hamiltonian has a simple transformation property under U- 
spin. Its three pieces responsible for Cabibbo-favored (CF), singly Cabibbo-suppressed (SCS) 
and doubly Cabibbo-suppressed (DCS) decays transform, when normalized suitably, like 
three components of a U-spin triplet operator denoted (U = 1, 1/3 = —1, 0, +1): 


ttCF _ 

n w — 


H, 


SCS 

w 


Hy 


DCS 

W 


G 

cos 2 9c(s c ) (ud) = — cos 2 0^(1, —1) , 

\/2 

Gp - 

—j= cos 9c sin 9c[(s c)(u s) — ( dc){ud )] = \/2cos9c sin0c(l, 0) , 
V 2 

G F , 2 
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sin 2 9c(dc)(us) = — sin" 9 c (l, +1) . 


(5) 


We have suppressed the chiral structure of V-A operators, using V u d = V cs = cos 9c,V us = 
—V c d = sin 9c for Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. 

Charmed meson decays are dominated by real CKM matrix elements associated with the 
first two quark families. For the most part we will be using a parametrization of the CKM 
matrix up to terms which are fourth order in A = V us HU. In Section ED we will discuss the 
effect of higher order terms in A on four ratios of amplitudes studied in this section. 

Virtual b quarks in penguin amplitudes may lead to tiny direct CP asymmetries in SCS 
decays of order (a s (m 2 c )/p){\V c bV u b\/\V cs V us \) ~ 10~ 4 p2], depending on the final state. No 
CP asymmetries of this order are expected in CF and DCS decays. An order of magnitude 
larger asymmetries may occur in SCS decays which are subject to potential penguin ampli¬ 
tude enhancement [19]. No CP asymmetries at these small levels have been measured so 
far. We will seek amplitude relations which hold at this high precision but not at higher 
accuracy. Errors in amplitudes are half the errors measured in decay rates. Thus we will 
neglect in our discussion direct CP asymmetries, assuming that branching ratios for D° —y f 
are given by averages measured for decay processes and their CP conjugates, 


B(D° -+f) = B(D° /) CPav = -+f) + B(D° /)] . (6) 


The D° is a U-spin singlet. Denoting hadronic matrix elements in the U-spin symmetry 
limit by a superscript (0) and defining a reduced matrix element, A = (1, t/ 3 | (1,1/ 3 )|0, 0), 
one obtains from Es. (j3j) and ([5]) 


{n+K-\H° F \D°)W 


cos 2 9q 


A , 


(7) 
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(K+K- - n + 7r-\H^ s \D 0 )W 
cos 9 C sin 9 C 
(K+7t-\H^ cs \D 0 )^ 


= 2 A 
= A . 


— sin 9 C 

Eq. (HJ) leads to 

(K+K~ + n + TT-\H^ s \D°) {0) oc (0, 0|(1, 0)|0, 0) = 0 . 
Using © this implies 

(K+K~ | H$f s | D °) (0) (tt+tt- | H^ s | D °) (0) 


cos $c sin 9 C 


— cos 9c sin 9c 


= A . 


( 8 ) 

( 9 ) 


( 10 ) 


( 11 ) 


Thus the four amplitudes in (J7J), (j9]) and (HIT) denoted by the decay final state, A(f) = 
(f\Hw\D°) have simple ratios in the U-spin symmetry limit [I], 


A°\'k+K ) : A^\K+K ) : A <K °\'K+'n ) : A^\K+ti ) = 1 : tan 9c ■ — tan^p : — tan 2 . 

(12) 

We note that a derivation of the two ratios, A°\K+it )/A(°\it+K ) = — tan 2 dc and 
A°\7t+ti-)/A^(K+K-) = -1, uses only symmetry under d <(->■ s reflection ra implying, 


(ir+K |(sc)(ud)|-D 0 ) = (Ii+ir \(dc)(us)\D°) , 

(. K+K~\(sc)(us ) — ( dc)(ud)\D °) = — (7r + 7r _ |(sc)('us) — (dc)(ud)\D°) , (13) 


and does not require full SU(2) U-spin symmetry. This full symmetry is required for the 
relation between these two pairs of processes. 


IV First order U-spin breaking 

First order U-spin breaking corrections to the amplitudes ((7j), (j9j) and OB are obtained 
by multiplying the weak Hamiltonian or the final states by an s—d spurion mass operator, 
Mubrk oc (ss) - (dd) [2j: 

(f\H eS \DY ] = (f\H eS M Vhlk \D°) + ( M Vhrk f\H eS \D°) . (14) 


While for CF and DCS decays one has simply 


rrCF,DCS7i f _ rrur J «■ 

n eS 1V1 Ubrk — n w 1V1 Ubrk 


(15) 


the effective Hamiltonian for SCS decays obtains at first order an additional nonperturbative 
s + d penguin term P s+ d due to an s — d mass difference 


Ffeff S ^Ubrk — H^ S Mubrk + P S +d 


l w 


(16) 
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The U — 0 penguin amplitude in SCS decays interferes with opposite signs with the U — 1 
tree amplitudes in D° —> K + K~ and in D° —> 7i + 7r~. This may potentially increase the first 
amplitude and decrease the second. This effect of the penguin amplitude has been pointed 
out very early in Refs. [2T[ [22, 23], and has been studied recently in Refs. [20, |2H 125] with 
its implication on CP asymmetries in these processes. 

Eqs. (I15p and (IIBI) imply different first order U-spin breaking corrections in D° —> tt + K~ 
and D° —>■ K + tt~, on the one hand, and in D° —> K + K~ and D° —> on the other. 

We will show now that the corrections within each of these two pairs of processes have equal 
magnitudes and opposite signs when normalized by the corresponding U-spin symmetric 
amplitudes. We will first follow a full SU(2) U-spin argument presented in [9], and then 
derive this result in a simpler manner using d t-?- s reflection. 

Consider first Eq. (ITTjl for decays to the two U-spin triplet states, |/i) = \i t + K~) = 
— 11, —1) and |/ 2 ) = \K + n~) = 11,1), to which (TT5]l applies. and H^p s transform like 
(1, —1) and (1, +1) while the s—d spurion mass operator Mubrk behaves like (1, 0). Since the 
D° is a U-spin singlet only the triplet operators in the products ff^ CS,CF Mubrk oc (1, ±1)(1, 0) 
contribute to the triplet final states ±|1, ±1), and only the triplet states in M Vbrk \l, ±1) oc 
(1, 0) 11, ±1) obtain contributions from the triplet Hamiltonian operator. Thus the two terms 
in (ITi]) involve coupling of the product (1, ±1) <g) (1, 0) into (1, ±1), where ± signs correspond 
to K + tt~ and n + K~ final states. Using a property of Clebsch-Gordan coefficients, 

(l,l;n, 0|1,1) = (—l) n (l,—l;n, 0|1,—1) , 

we therefore find 

(n + K-\H^\D 0 )M _ (K + n-\H^ cs \D 0 }^ 
cos 2 6c sin 2 9q 

That is, first order U-spin breaking terms contribute equally but with opposite signs to D° —y 
7 t + K~ and D° —> K + n~, when normalized by corresponding U-spin symmetric amplitudes 
given in Eqs. ([7]) and ([9]): 

{T+K-\HS, F \Dy» (A-+>r-|g° cs |r>°)W _ 

We denote by e\ the U-spin breaking term in D° —> K + n~ normalized by its U-spin invariant 
amplitude. 

A simpler and rather immediate derivation of (TT9]l may be obtained by applying d -h- s 
reflection to Eq. (1141) for / = tt + K~. Noting that the s — d spurion mass operator Mubrk 
changes sign under this reflection, one has 

(7r + W-|(sc)(hd)M Ubr k| J D 0 ) + (M Vhrk (n + K-)\(sc)(ud)\D°) 

= -(K + n-\(dc)(us)M Vhlk \D 0 ) - (M Vhrk (K + n-)\(dc)(us)\D 0 ) . (20) 

This leads directly to (flSll and (flOl) . 

This short argument applies also to first order U-spin breaking in SCS decays, D° —y 
K + K~ and D° —> 7r + 7r _ , since the penguin operator P s +d in (TT6l) also changes sign under 


(17) 

(18) 
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d -H’ s reflection. Therefore, 


fr + *-|*& c3 |P°) (1) (A-+A-|g|C s |£>°)W _ 

(ir+!T-|i/g, cs |B»)(») (Jl+/l-|ffg? s |D»)(») 2 ’ { ’ 

where e 2 denotes the U-spin breaking term in D° —y K + K~ normalized by its U-spin invariant 
amplitude. Furthermore, the change in sign between first order terms in D° —> K + K~ and 
D —>■ 7 r + 7 r _ applies separately to contributions of the two operators on the right-hand side 
of (fT 6 i) representing tree and penguin amplitudes. This leads one to expect that the U-spin 
breaking parameter 62 , involving both tree and penguin amplitudes, is larger than ei which 
involves only tree amplitudes [ 20 j. 

Combining the results (TTUf) and (OTj) with the zeroth order Eqs. (J7J), © and (TTT1) one 
obtains the following first order expressions for decay amplitudes: 


A(D° 7 t + K~) 

A(D° ->• K+ tt“) 
A(D° —> 7i + 7r ~) 
A(D° -> K + K~) 


cos 2 9 C A(1 — ei) , 

— sin 2 9c A{1 + ei) , 

— cos 6c sin 6 C A( 1 — e 2 ) , 
cos 9c sin 9c A( 1 + e 2 ) • 


( 22 ) 


V U-spin breaking of arbitrary order 


U-spin breaking of order n in decay amplitudes (f\H c g\D°) is obtained by introducing in 
the Hamiltonian or in the final state a total of n powers of the s — d spurion mass operator, 
applying (ITTH) to SCS decays. Generalizing the argument for a relative negative sign in first 
order breaking, based on a change of sign of Mubrk and P s +d under d t-?- s reflection, we 
conclude that a negative relative sign applies to odd n and a positive sign to even n: 


(tt + K-\H^\D 0 )^ 
cos 2 9 C 

(K + K-\H$f s \D 0 )^ 
cos 9c sin 9c 

Thus we may sum up: 


(-ir 

= (-1) 


(K + n-\H^ s cs \D 0 )^ 
— sin 2 9c 

n (lT + TT-\H$f S \D 0 )^ 
— cos 9c sin 9c 


A(D° ->• 7 t + K~) 
A(D° ->• K + n~) 

A(D° —> 7r + 7r _ ) 

A{D° K + K~) 


cos 2 9 c A[1 - e 1 + aie 2 - a[e\ + ...] , 

— sin 2 9c A[ 1 + + d\t\ + o!^e\ + ...] , 

— cos 9c sin 9c A[1 — e 2 + a 2 e \ — a' 2 e\ + ...] , 
cos 9c sin 9c A[~L T e 2 T o 2 e 2 -f- u 2 e 2 T ...] . 


(23) 

(24) 


(25) 


While the complex U-spin breaking parameters ei i2 and the nonerturbative coefficients 
°i, 2 , a i 21 ■■■ are n °t calculable from first principles, one expects the first two parameters to be 
around 0.2 — 0.3 (e 2 being larger in magnitude than for the above-mentioned argument) 
and the coefficients to be of order one, |ai i2 | ~ \a\ 2 | ~ 1. 
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Expanding magnitudes of amplitudes up to and including third order we find 

|1 ± e + ae 2 ± a'e 3 \ = 1 ± Ree + -(line ) 2 + Re(ae 2 ) ± Re(a'e 3 ) 

=F -Re e (Im e ) 2 ± Im e Im(ae 2 ) . (26) 

In the next section we will use this expansion for studying ratios of magnitudes of decay 
amplitudes, identifying relations among ratios in which U-spin breaking terms up to third 
order cancel, thus being sensitive to tiny fourth order U-spin breaking. We will then argue 
that there is no need to go beyond third order in (|26li for showing that fourth order terms 
do not cancel in these relations. 


VI High-precision relation among ratios of amplitudes 


We define four ratios of amplitudes: 


Ri 

R 2 

r 3 

i?4 


I A(D° ->■ K+ 7T-)| 

\A{D° 7r + A' _ )| tan 2 6 C ’ 

| A(D° ->■ K+K ~)| 

\A(D° ^7T+7T-)| ’ 

\A(D° ->• K+K~)\ + \A(D° -> 7T+7r-)| 

| A(D° —> 7t + A'~)| tan 9c + \A(D G —> /v+7r _ )| tan -1 6 c ’ 


N 


\A(D° -»• K+K~)\\A(D° ->■ 7r+7r-)| 
\A(D° 7t + K~)\\A(D° -> K+ tt-)| ' 


These four ratios are not mutually independent. They obey a trivial identity, 


f?4 




l - [(fi 2 


\ i - [(ft 


i)/(ft + ig 

i)/(ft + i ) 2 ' 


(27) 


(28) 


We will prove a nontrivial relation involving R 3 and R 4 and another nonlinear function of 
Ri and R 2 , that holds up to maximal fourth order U-spin breaking. 

Using (126)1 we start by expanding the two ratios R\ and i? 2 U P to third order U-spin 
breaking: 


R\ — 1 + 2[Re 64 + (Re ei) 2 ] + (9(e 3 ) , 

R .'2 = 1 + 2[Re £2 + (Re e 2 ) 2 ] + 0(e 2 ) • (29) 

These two ratios involve first order corrections given by 2Reei and 2Ree 2 . Second order 
corrections in these ratios are given by squares of these same real parts with no dependence 
on the unknown parameters Gp and a 2 . Thus measurements of R\ and R 2 provide a way for 
calculating Re E\ and Re e 2 up to third order corrections. Solutions for Re and Re e 2 using 











Eqs. (125]) should include the U-spin symmetry limit, requiring Reei = 0 and Ree 2 = 0 for 
Ri — 1 and i ? 2 = 1, respectively, rather than Reei = —1 and Ree 2 = —1. This implies 


Reei — ^ ~ + ^( e i) > 

Ree 2 = - (^]/2R 2 — 1 — 1^ + 0(e%) . (30) 

As we will see immediately, these two first order U-spin breaking parameters do not determine 
only Ri and i ? 2 but also the difference R 3 — R 4 . 

The two ratios R 3 and R 4 , in which Erst and third order terms cancel [26l 127] , may be 
expanded up to fourth order: 

R 3 = 1 + 7 ^[(hn e 2 ) 2 - (Imei) 2 ] + Re [a 2 e\ - cpe 2 ] + 0(ef,e%) , 

Ra = 1 — -Re (e 2 — e 2 ) + Re (o 2 e 2 — Gqe 2 ) + C?(e 1 , 62 ) 

= 1 + ^[(Ime 2 ) 2 - (Imq) 2 ] + Re (a 2 e 2 - eqe 2 ) - ^[(Ree 2 ) 2 - (Reei) 2 ] 

+ o(4,4)- (si) 


These two ratio are noticed to differ by second order U-spin breaking terms depending solely 
on Re ex and Ree 2 : 

R 3 — Ra = -[(Ree 2 ) 2 — (Reei) 2 ] + 0(ef, e t) ■ (32) 

Using 

(Ree ,) 2 = i (j/2Ri - 1 - l) 2 +2Re Ci 0(e?) = ^ (^2 R { - 1 - +0(e$) , * = 1,2, (33) 

one obtains the the following nonlinear relation among the four ratios of amplitudes, which 
holds up to tiny fourth order U-spin breaking terms: 


AA = R3 — Ra -\— 



0(4,4 


(34) 


This relation may also be obtained by expanding the identity (128]) to first order in [(i?i )2 — 
1 )/(-Ri ,2 + 1 )] 2 and applying (]29]) and (]30]) . Fourth order terms in (Ree *) 2 are proportional to 
Reej. It can be easily shown that this is not the case for terms of this order occurring in A 3 , 
Ra and in their difference. Therefore one concludes that fourth order terms do not cancel in 
A R and in any higher order expansion of the square root in (]28j) . Thus, in hindsight, there 
is no way of obtaining a U-spin breaking relation of higher order than (152]) and one does not 
need to go beyond third order in (|26|) . 

In our derivation of (152]) we have used a parametrization of the CKM matrix up to terms 
which are fourth order in A = V us = 0.2253 ± 0.008 [T7]. We now study the effects of higher 
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order terms in A on Ri and on this relation. Including A 4 and A 5 terms in the CKM matrix 
one has [23] 


v ud 

= l - -a 2 - -A 4 , 

2 8 ’ 


Vus 

= A, 


V cd 

= -X + ^A 2 X 5 [l-2 (p + iri)], 


V cs 

= l - Ia 2 - -A 4 - -A 2 A 4 , 

O Q 0 ’ 

(35) 


where A = 0.82 ± 0.02, p = 0.12 ± 0.02, rj = 0.36 ± 0.02 [IT]. 

We define a small parameter £ determining the effect of A 4 and A 5 terms on ratios of 
amplitudes: 


i + £ = 


\A 2 A 4 |l - 2{p + irj)\ 
1-±A 2 A 4 


(36) 


Using the above values of A ,A,p and r) we find |£| ~ 10” 4 . While R 4 is unaffected by £, R 3 
and the two quadratic terms in ([34]) obtain corrections suppressed by £ and by first order 
U-spin breaking. These corrections, respectively — |£(Reei + Ree 2 ), + |£Reei and +|£Ree 2 , 
cancel each other. The remaining corrections in (I34[) . suppressed by £ and by second order 
U-spin breaking, are much below the level of 10~ 4 and may be safely neglected relative to 
tiny fourth order U-spin breaking terms of order 10~ 3 or 1CU 4 . 

Taking e, ~ 0.2 as a typical value for first order U-spin breaking, the nonlinear relation 
A R = 0 is expected to hold at a very high precision of order e 4 ~ 10~ 3 . We will confirm 
this prediction in Section IVIII1 At this high precision one cannot ignore first order isospin 
breaking terms, which one would generally assume to be around (m d — m u ) Aq CD ~ 10” 2 . 
We will study these terms in the next section showing that, in fact, they are suppressed by 
both isospin breaking parameters and by second order U-spin breaking. 


VII First order isospin breaking 


First order isospin breaking is introduced by multiplying the weak Hamiltonian by a d — u 
spurion mass operator, 


MjhrV oc 


(rlrl 




(37) 


transforming like a combination of a U-spin singlet and triplet. Isospin breaking contributions 
of the U-spin singlet operator in the four amplitudes (1251) are identical when normalized by 
suitable CKM factors, and may be absorbed into the U-spin symmetric amplitude A. This is 
true also for U = 0 isospin breaking electromagnetic interactions because the d and s quarks 
have identical charges. 

Contributions of the triplet operator in (137p follow the signs of first order U-spin breaking 
corrections. They are represented by two distinct first order isospin breaking parameters, 
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<5i - for U-spin triplet states n + K and K + tt , and h 2 - for K + K and 7 r + 7 r , the two 
component of a U-spin singlet state. Thus 


A(D° ->• 7 t + K~) 
A(D° -+ K + tt~) 
A(D° —> 7r + 7r _ ) 
A(D° -+ iT+AT-) 


cos' 6q A(l — £i T cqe 2 — cqe 3 — hi T ...) , 

— sin' 9c A{\ + Ci + cqe 2 + eqe 3 + hi + ...) , 

— cos0csin0c A(1 - e 2 + a 2 e 2 + —a 2 e 2 - h 2 + ...) , 
cos 9c sin 6c A{\ + e 2 + ® 2 C 2 H—bo 2 e 3 + h 2 + ...) . 


(38) 


Instead of (1261) we now expand: 


1 ± e + ae 2 ± a'e 3 ± 5 \ 


= 1 ± Re e + -(Ime ) 2 + Re(ae 2 ) ± Re(a , e 3 ) 

+ —Re e (Im e ) 2 ± Im e Im(ae 2 ) ± Re h + Im h Im e . 


(39) 


This characteristic amplitude expansion includes two new terms, ±Re h and +Im h Im e, the 
latter involving suppression by both isospin and U-spin breaking parameters. We will show 
that terms of this order do not affect the nonlinear relation (I34j) . 

The expansion of the four ratios of amplitudes now includes terms which are first order 
in isospin breaking and other terms suppressed by both isospin and U-spin breaking: 


Ri = 1 + 2[Ree! + (Reei) 2 ] + 2Rehi + 4Rehi Reei + (9(e 3 ) + (9(hie 2 ) , 

R 2 = 1 + 2[Re e 2 + (Ree 2 ) 2 ] + 2Reh 2 + 4Re h 2 Re e 2 + 0(e 2 ) + (9(h 2 e 2 ) , (40) 


#3 = 1 + -[(Ime 2 ) 2 - (Imei) 2 ] + Re (a 2 e 2 - cpe 2 ) + (Im S 2 Ime 2 - Erich Imei) 

+ ^(hlg) + ^(hl,2Cl ;2 ) ) 

#4 = 1 + -[(Ime 2 ) 2 — (Imei) 2 ] + Re (a 2 e 2 — eqe 2 ) + (Imh 2 Ime 2 — Im hj Imei) 

— — [(Re e 2 ) 2 — (Re ei) 2 ] — (Re h 2 Re e 2 — RehiRe ei) + (A(ef 2 ) + ^(hige 3 ^) • (41) 


Eqs. (1401) imply for i — 1, 2 
1 


Ree, = 


(^\j2Ri — 1 — 1 ^ — Rehj — 2 RehjReej + 0 (6^) + 0 (e 3 ) , 


(42) 


or 


(Ree 8 ) 2 = - (yj2Ri - 1 - lj - 2Reh^Ree 8 + 0(84) + O(ef) . (43) 

Eqs. (j4Tj) lead to 

#3 — #4 = — (Re e 2 ) 2 — Reei ) 2 + (Reh 2 Ree 2 — RehiReei) + 0(ef 2 ) + ^(hi, 2 e 3 2 ) • (44) 
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Consequently isospin breaking terms of the form Re5 2 Ree2 — Re^Reei cancel in A R: 


A R = 


R,—Ri 



= 0{e\ ) e A 2 )+0{8 l el5 2 el). (45) 


It is remarkable that isospin breaking modifies the nonlinear relation (154)1 by terms which 
are suppressed by both first order isospin breaking and second order U-spin breaking factors. 
Taking Si ~ 10 -2 , e* ~ 0.2 — 0.3, these terms are expected to be at most 10 -3 , similar in 
magnitude to fourth order U-spin breaking terms affecting this relation. 


VIII Experimental tests in D° — y P + P 

We will now apply current experimental data to study the hierarchy among U-spin breaking 
terms of increasing order. Our final goal is testing the predicted amplitude relation (145)1 . 
Hadronic decay amplitudes (A) are obtained from measured branching ratios (J3) by elimi¬ 
nating phase space factors depending on final particles center-of-mass 3-momenta (p), and 
on the D meson mass and its lifetime (Mp and Tp), 

\A\ = M D *[^E . (46) 

V t dP 

In our calculation of amplitudes we will disregard a common factor Md^JPkJtd which cancels 
in the four ratios Ri. Values for measured branching ratios ra, center-of-mass momenta, 
and amplitudes defined in this manner are quoted in Tables HI 

Note that all four amplitudes include a factor corresponding to the branching fraction 
of the Cabibbo-favored decay D° — > ir + K~. We have included no error in the amplitude 
for this process as the other three branching ratios (including errors) have been measured 
relative to this process mt The relative errors in the amplitudes of these three processes 
are all below the level of one percent. We will assume no correlation between these errors, 
which have been measured in three independent analyses for different final states. The high 
precision achieved recently by the CDF, LHCb and Belle collaborations in measuring the 
DCS amplitude is remarkable [29], as it required time-dependent separation between this 
highly suppressed decay and D°-D° mixing followed by the CF decay. 

Using values of amplitudes given in Tables Q] and tan 9 c = 0.23125 ± 0.00082 [17] we 
calculate the four ratios Ri defined in Eq. (127)) . 

R l = 1.115 ±0.012, 

R 2 = 1.811 ±0.020, 

R 3 = 1.052 ±0.008, 

R 4 = 1.008 ±0.007. (47) 

Errors in the ratios have been obtained by adding in quadrature errors in the relevant 
amplitudes. 
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Table I: Branching fractions and amplitudes for D° —» P + P decays [IT] 


Decay mode 

Branching fraction 

p (GeV/c) 

\A\ = ^jB/p (GeV/c) -1//2 

D° 

7 t+K- 

B wK = (3.88 ± 0.05) xl0“ 2 

0.861 

1.078B^ 2 

D° 

K+ 7T“ 

(3.56±0.06)x10- 3 H^ 

0.861 

(0.06430 ± 0.00054)H^ 2 

D° 

—> 7T + 7r _ 

(3.59 ± 0.06) xlO" 2 ^ 

0.922 

(0.1973 ±0.0016)H^ 2 

D° 

-> K+K~ 

(10.10 ±0.16) xlO- 2 H^ 

0.791 

(0.3573 ± 0.0028)H^ 2 


We have seen that the first two ratios involve first order U-spin breaking terms. These 
terms, depending on two distinct U-spin breaking parameters e 4 and 62 , are considerably 
larger in R 2 than in R\. This has been anticipated in the discussion below Eq. (12 1 p . Using 
(02]), in which we neglect first order isospin breaking and third order U-spin breaking, we 
calculate reasonably small U-spin breaking parameters, 

Reei = 0.054 ±0.005, 

Re e 2 = 0.310 ±0.006. (48) 

The other two ratios, R 3 and R 4 , given in (011) in terms of ei )2 and coefficients gr 2 of order 
one, deviate from one by second order U-spin breaking terms (we neglect terms suppressed 
by both isospin and U-spin breaking and fourth order terms in U-spin breaking), 

R 3 — 1 = ^[(Ime 2 ) 2 — (Irnei) 2 ] ± Re {a 2 e\ — Gpe^) = 0.052 ± 0.08 , 

R 4 - 1 = Re[(a 2 - ^)e 2 2 - (aj - ^)ef] = 0.008 ±0.007. (49) 

The hierarchy between the first order parameters in and the second order terms calculated 
in [jDj confirms and justifies the perturbative approach we have applied in this study to U- 
spin breaking. Without having a method for calculating the nonperturbative coefficients a,, 
the almost exact cancellation of second order terms in i ? 4 seems to be accidental. In view 
of the small value of Re ei and the much larger value of Re e 2 this approximate cancellation 
seems to imply a 2 ~ 1 / 2 . 

Having shown that second order U-spin breaking terms are a few percent, one expects 
fourth order terms to be of order 10 -3 , comparable in magnitude or larger than terms which 
are first order in isospin breaking and second order in U-spin breaking. Let us now check 
this prediction in the relation (145|) which contains on the right-hand side terms of these two 
kinds. Using the definitions of R t in (1771) and adding in quadrature errors in amplitudes 
[rather than errors in Ri given in (1471) ]. we obtain 

AR = (-3.2 ± 0.4) x 10 -3 . (50) 


This confirms our prediction. 


13 








IX Experimental tests in D° —> V + P 


As mentioned, the discussion in Sections HHI - IVIII applies also to three other classes of 
processes involving one or two charged vector mesons, D° —> V + P~,D° —> P + V~ and 
D° —> V + V~. In particular, a nonlinear amplitude relation similar to (|45l) holds in each one 
of these classes with a precision depending on the size of U-spin breaking. In this section 
we summarize concisely the situation relevant to this question in D° —» V + P~, consisting of 
the four processes, D° —> p + K~ , D° —>■ K* + n~, D° —> p + 7T~ and D° —> K* + K~. We denote 
first order U-spin breaking and isospin breaking parameters in these processes by e[ 2 and 
5[ 2 , respectively, in analogy to e lj2 and S 1:2 in D° P + P~. 

Dehning four ratios of amplitudes R[ in analogy with (l27jh 


R[ 

R-2 

R's 


R 


/ 

4 


| A(D° -> AT*+7t-)| 

\A(D° —>■ p + K~) \ tan 2 6 C ’ 

| A(D° ->■ K*+K~) | 

\A(D° —)• p + 7i~)\ ’ 

| A(D° -)■ K*+K~)\ + | A(D° ->■ p+7T-)| 

|A(D° ->■ p+/l~)| tan0 c + l^(^° ->• K* + tt~)\ tan" 1 ’ 


\ 


\A(D° ->■ Ji | |A(U>° ->■ p+7T-)| 
|A(D 0 ->■ p+A'-)||Ap° ->■ JW+7T-)| 


5 


(51) 


one obtains a sum rule analogous to (l45|) : 


A R’ = R’-R’+- 


2-R'i -1-1 - -1-1 


= 0(e' 1 4 ,e' 4 ) + 0(5' 1 e^) 


(52) 

Magnitudes of amplitudes for the P-wave decays D° —> V + P are obtained from corre¬ 
sponding branching ratios using 


\A\ = M D \ 


I 8ttI3 
tdP 3 


(53) 


We will disregard again the factor M d ^8it/td since we are only concerned with ratios of 
amplitudes. Values for measured branching ratios mi. center-of-mass momenta, and am¬ 
plitudes defined in this manner are quoted in Tables HU Relative errors in CF and SCS 
amplitudes are reasonably small, between two and three percent. In contrast, the relative 
error in the DCS amplitude | A(D a —> K* + 7 r _ )|, obtained by the CLEO and BaBar collabo¬ 
rations through Dalitz plot analyses of D° —> Ks7 r + 7r [30], is quite large, This large 

asymmetric error limits considerably the precision of R[, R' 3 and R' A . Adding in quadrature 
errors in relevant amplitudes, measured independently for different three-body final states, 
we calculate 


K = 


o.97i±8:2£ 
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Table II: Branching fractions and amplitudes for D° —y V + P decays pTJ 


Decay mode Branching fraction p (GeV/c) \A\ = JB/p 3 (GeV/c) 3 / 2 


D° p+K~ 

0.108 ±0.007 

0.675 

0.593 ±0.019 

D° -)■ K*+ 7T" 

(3.42t]; 3 °)xl0- 4 

0.711 

0.0308™ 4 

D° —>■ p + 7T~ 

(9.8 ± 0.4) xlO -3 

0.764 

0.148 ±0.003 

D° K*+K- 

(4.38 ± 0.21)xl0” 3 

0.610 

0.1389 ±0.0033 


R' 2 = 0.939 ± 0.029 , 

R' 3 = 1.061™, 

K = 1-0611^42 ■ (54) 

Expansions similar to (TT0]1 and (PTTj) apply to these four ratios in terms of e' x 2 and 2 . 
The leading U-spin breaking corrections in R[ and R 2 are hrst order, while in R 3 and R’ A 
they are second order. The measured values of the first two ratios imply 

Re e', = j - 1 - l) = -0.015™ , 

Ree'j = 1 Uwi - 1 - li = -0.032 ± 0.016 . (55) 

We note that the numerical value of Ree 2 is smaller by about an order of magnitude than 
the value of Re 62 calculated in (l48lh That is, first order U-spin breaking in SCS D° —$■ 
V + P~ decays is about an order of magnitude smaller than in corresponding D° —» P + P~ 
decays. This seems to imply that an enhancement of the U-spin breaking penguin amplitude 
suggested to occur in D° —> P + P~ [19] is not at work in D° —>■ V + P~. The other U-spin 
breaking parameter, Ree' 1; does not involve a penguin amplitude. In spite of the current 
large error in this parameter the first of Eqs. (157)]) favors strongly iRee'J < 0.1, suggesting 
that a value close to that measured for |Reei| is not unlikely. 

These values of Ree' 12 imply that typical second order U-spin breaking terms in R' 3A 
are around one percent. Confirming this prediction and obtaining a more precise value for 
Ree^ requires a substantial improvement in the measurement of B(D° —>■ K* + n~). I11 the 
meantime we use the values measured for iRee^ 2 | to argue that fourth order U-spin breaking 
terms should be around 10~ 4 . This is also expected to be the magnitude of terms in A R' 
suppressed by both isospin breaking and by second order U-spin breaking. Thus, due to 
smaller U-spin breaking parameters in D° —>■ V + P~ relative to D° —>■ P + P~ one predits 
A R' to be about an order of magnitude smaller than A R. Using the amplitudes in Table [TT| 
we calculate 

A Ft = (0.2+ 3 ; 2 ) x 10" 4 , (56) 

where the error is dominated by the error in \ A(D° —> K* + tt~)\. This confirms our prediction. 
A more precise test could be achieved by improving the measurement of the branching ratio 
for this process. 
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Table III: Branching fractions and amplitudes for D° —y P + V decays [TTj 


Decay mode 

Branching fraction 

p (GeV/c) 

\A\ = ^B/p 3 (GeV/c)" 3 / 2 

D° ->■ 7 t+K*~ 

(4.98l^i) xlO" 2 

0.711 

H 070+U.UI7 

u.cwz_ 0019 

D° K + p~ 

- 

0.675 

- 

D° —> TT + p~ 

(4.96 ±0.24)xlO” 3 

0.764 

0.1055 ±0.0026 

D° Jl+Jl*- 

(1.56 ± 0.12)xlO” 3 

0.610 

0.0829 ±0.0032 


X Experimental tests in D° —> P + V 


Current branching ratios and amplitudes for D° —>■ P + V~, consisting of D° tt + K*~ , K + p ~, 
7r + p _ , K + K*~ are given in Table III. Ref. [T7] quotes no branching ratio measurement for the 
DCS decay D° —> K + p~. Three-body decays, D° —> /i + 7r _ 7r°, involving K + p~, K* + n~ and 
other intermediate states, have been observed by Belle [31] and BaBar [32] with branching 
ratios (3.18 ±0.29) x 10 -4 and (2.97±0.19) x 10” 4 , respectively. The D° —> K + 7r~7r° events 
involve interference of DCS decays with D°-D° mixing followed by CF decays. Evidence 
for D°-D° mixing at 3.2 standard deviation was presented by Babar [33], measuring the 
fraction of K + p~ in these events to be (39.8 ± 6.5)%. (No interference would have implied 
B(D° —>■ K + p~) ~ 1.2 x 10~ 4 .) More work is needed for resolving the effect of D°-D° mixing 
on these events, and for obtaining a solid measurement of B(D° —> K + p~). 

We define ratios of amplitudes R" in analogy with (137]) . 


K 

K 

K 

K 


\A(D° -)• K+p~)\ 

\A(D° —>• n + K*~)\ tan 2 9 C ’ 

\A(D° ->■ K+K *~)| 

\A(D° —> 7T+p _ )| 

\A(D° -> K+K*~) | ± | A(D° -> 7r+p-)| 

| A(D° tt+K*-)\ timOc ± \A(D° K+p~)\ tan' 1 9 C ’ 


N 


-> K+I<-)\\A(D a -> ir+p~)l 
\A(D° ->jr+A'-)P(£>» -> K+fr )I ’ 


(57) 


We denote first order U-spin breaking and isospin breaking parameters in these amplitudes 
by e" 2 and S'{ 2 , respectively. 

In the absence of a solid measurement of |^4(Z1 0 —> K + p~)\ one can only calculate R 2 . 
Neglecting isospin breaking and third order U-spin breaking one obtains 


R 2 = 1 ± 2 


Reed 


± (Re e 


n\2 


= 0.786 ± 0.036 


(58) 


which implies 


Reed = -0.122 ±0.024 . 


(59) 
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That is, the magnitude of the U-spin breaking parameter in SCS D° —> P + V~ decays is 
intermediary between corresponding parameters in D° —> P + P~ (Re 62 = 0.310 ±0.006) and 
D° —)■ V + P~ (Ree' 2 = —0.032 ± 0.016). Namely, no significant U-spin breaking penguin 
enhancement applies to D° —>■ P + V~. One expects a similar or smaller magnitude for Ree". 

The two ratios R" and R " deviate from one by second order U-spin breaking terms [see 
Eqs. (El])] which are expected to be at most a few percent. Using R% — 1 ± 0.05, where we 
include a conservative uncertainty of 5% due to second order U-spin breaking corrections, 
one obtains the following prediction for B(D° —> A' + p _ ) (34] . 

\A(D° —► K + p~)\ = 0.0237±0.0025 => B(D° K + p~) = (1.73 ± 0.36) x 1(T 4 . (60) 

For comparison, assuming R" — 1 ± 0.05 implies a very similar prediction, 

| A(D° ->• K + p~) | = 0.0235 ± 0.0028 => B(D° -> K+p~) = (1.70 ± 0.40) x 1(T 4 . (61) 


This value of \A(D° —> A' + p _ )| would imply Ree'/ = 0.08 ± 0.06, comparable in magnitude 
to Re e '2 and in agreement with expectation. 

Taking for \A(D° —> K + p~)\ the value in d6T]) . using the three measured amplitudes 
quoted in Table III, and assuming no error correlation between the four amplitudes, one 
obtains 


A R" = R'n - A" ± - 

S 4 g 



= (0.8l“) x IQ' 4 . 


(62) 


This value, which is similar to (15^1) . is in agreement with our prediction that, in view of 
the above values of Ree'/ 2 , fourth order U-spin breaking terms and isospin breaking terms 
suppressed by second order U-spin breaking should be of order 1CT 4 . The positive error in 
A R" is dominated by the uncertainty assumed in the unmeasured DCS amplitude. One still 
awaits a solid measurement of B(D° —> A' + p - ) which would test the prediction (16TT) . The 
larger negative error in A A", originating in errors on the three measured amplitudes in Table 
III, may be reduced by improving the corresponding branching ratio measurements. 


XI D° decays to pairs of neutral pseudoscalars 

In the U-spin symmetry limit one obtains simple amplitude relations for D° decays to pairs 
of light neutral pseudoscalar mesons [35]. We will go through the symmetry argument first, 
extending it to include U-spin breaking at arbitrary order in CF and DCS decays. We will 
then demonstrate a few amplitude relations which hold up to second order U-spin breaking. 
Other relations of this kind have been studied in Ref. [8J in the framework of flavor SU(3). 

I 11 the symmetry limit one neglects p— 7 / mixing which is due to first order U-spin breaking 
represented by the spurion mass operator Mubrk- Thus, we will write amplitudes for r] = p 8 
in our discussion of the symmetry limit, while r/ 8 will be used explicitly when introducing 
U-spin breaking. When discussing amplitudes, rates and asymmetries for p 8 we will assume 
(as has been assumed in Ref. [8j) knowledge of amplitudes including a relative phase for 7 ] 
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and rf and favored values of the mixing angle. Taking r/ = r/ 8 = (2ss — uu — dd)/y/ 6, the 
following superpositions of single neutral particle states belong to a U-spin triplet, 


| K°) = | ds) = |1,1) , \K°) = | sd) = -|1, -1) , 

while the orthogonal U-spin singlet is 


ky/3\v) ~ k°» = \ss - dd)/V2 = 11,0) , 

(63) 


-(|r]) + \/3|7r 0 )) = |ss + dd — 2uu)/V6 = |0, 0) . 


(64) 


Two-particle states in an S-wave are obtained by symmetrizing products of single particle 
states. 

Symmetrized products of two U = 1 single-particle states consist of U = 0 and U = 2 
(with U 3 = ±1, 0) states, while the product of two U = 0 states is pure U — 0. D° decay 
matrix elements of the U — 1 weak Hamiltonian (J5]) vanish for each one of these five states. 
This implies the following five U-spin symmetry relations: (For short notation we denote 
amplitudes by their final states.) 


V3A (0 \K 0 r]) -H (0) (A'°7r°) = 0, 

VSA i0 \K°r]) - A (0 \K°n 0 ) = 0, 

a(°\k°k°) = 0 , 

A^\rir]) + A^\ vrV) = 0, 

V3A w (t)tt 0 ) + V2A (0 \h°tt°) = 0. (65) 

Hadronic matrix elements for symmetrized products of U — 1, U 3 — ±1, 0 and U — 0 states 
are then given by a single {7 = 1 amplitude A: 

V3A (0) (K°rj) = A w (K°n 0 ) = -sin 2 9 C A , (66) 

V3A (0) (K°t]) = A w (K°tt°) = cos 2 6 C A, (67) 

VsA^' 1 (r]7T°) = V2 A^°\rjrj) = —V2A(°\ 7r°7r°) = \pi cos 6c sin 6c A . (68) 


Note that since we symmetrized final states also for identical particles, corresponding am¬ 
plitudes have been divided by \/2 in order that their squares give decay rates. 

The above amplitude expressions are analogous to the U-spin symmetry expressions ([7j), 
()9l) and dTTj) for decays into pairs of charged pseodoscalar mesons. Some of these relations 
follow merely from d -H- s symmetry. Since the {7 = 1 and {7 = 0 superpositions of 7 r° and ij 
states in (l63li and (l64]i are respectively antisymmetric and symmetric with respect to d f->- s 
reflection one has 


(K°(V3r) — 7i°)\(sc)(ud)\D°) = —(K 0 (V3r] — n 0 )\(dc)(us)\D 0 ), 

(. K°(ri + V37T 0 )\(sc)(ud)\D °) = (K°( V + V3n°)\(dc)(us)\D°) , (69) 
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implying 


a/3A(°) {K°rj) - A {0 \K 0 ^) 
cos 2 9 C 

A(°\K°ri) + V3 A(°\K°tt°) 
cos 2 9 C 


y/3A<®(K°ri) - A^(K°tt°) _ 

sin 2 9 C ’ 

AW(K°r]) + V3A( 0 \K 0 n 0 ) 4 

^9~c “73 


The right-hand sides follow using (l66j) and ()6T|) based on the full SU(2) structure of U-spin. 

U-spin breaking of order n in CF and DCS amplitudes is introduced by multiplying 
transition operators or final states by a total of n powers of the s — d spurion mass operator 
dTubrk which changes sign under d yy s. Consequently one has 


V3A( n \K°i l8 ) - A^\K°tt°) 
cos 2 9 c 

A^(K°r] 8 ) + y/3A^(K 0 n°) 

cos 2 9c 


y/3AW (It %) - AW 
sin 2 9c 

A^\K°r] 8 ) + y/3A^(K°n°) 
— sin 2 9c 


(71) 


Denoting first order U-spin breaking parameters in these two pairs of processes by eo and e' 0 , 
one may expand the above four linear combinations of amplitudes to arbitrary order, 


V3A(K°t] 8 ) - A(K°tt°) 
V3A(K°r) 8 ) - A(K°n°) 

^[A(K° m ) + V3A(K°n 0 )} 
-j-[A(K°ri 8 ) + V3A(K°n°)\ 


cos 2 9 c A[e 0 - a 0 el + ...] , 

— sin" 9cA[eo + + ...] , 

COS 2 0c^4.[l — 6g + UoCq 2 + •••] , 

— sin 2 9cA[1 + 6q + ciqc'q + ...] , 


(72) 


(73) 


where |a 0 | ~, |a()| ~ 1. 

Eqs. (I72p imply a linear amplitude relation in which first order U-spin breaking terms 
cancel, 

[V3A(K° V8 ) - A{K°tt 0 } tan 2 9 C + 73 A{K°r] 8 ) - A(K°n°) = 0 . (74) 

This relation has also been obtained using a general first order SU(3) breaking expansion |8], 
in which a dozen SU(3) breaking parameters contributing to these processes cancel in this 
relation. We note that, while it follows from Eqs. (1731) that the linear relation 

[A(K°r] 8 ) + VSA(K 0 7i 0 )\ tan 2 9 C - [. A{K°rj 8 ) + V3A(K { V)] = —= sin 2 9 C A (75) 

v 3 

is also free of first order U-spin breaking, the U-spin invariant amplitude A on the rignt- 
hand side is not necessarily invariant under the full flavor SU(3) group when including other 
states. This is just like the amplitude A defined above (J7J) for decays to charged particles. 

The four amplitudes for two body decays involving K° or K° and 7r° or rj 8 have the 
following first order expansions: 


A(K°tt°) 


cos 2 6c A 


-e 0 + 0(e, 


J2\ 


? 
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V^A(K°T]s) — COS' 9cA 1 — 6q + —6q + 0(Eq, £q 2 ) 


A(K 0 7t°) = — sin 2 6 C A 
V3A(K°rjs) = — sin 2 9c A 


1 + e o - 2 e ° + °( e o, e o ) 
1 + e o + ^ e 0 + e o) 


( 76 ) 


This implies that the two ratios of amplitudes in (1661) and fl67lh which are equal in the U-spin 
symmetry limit, are also equal when including first order U-spin breaking: 


Vs\A(K° V8 \ Vs\A(K° V8 )\ 


|A(A' 0 7t°)| 

\A(K°n°)\ 


\A(K°tt°)\ 

\A(K°r, 8 )\ 


— 1 + Re e 0 + 0[e 0 , e(,) , 


— 1 + 2Re e'o + 0(e 0 , e' 0 ) . 


(77) 


(78) 


| A(K ( °7T 0 ) | tan 2 6c \A(K°r] 8 ) | tan 2 6c 

Branching ratio measurements of the above two DCS decay modes are not feasible because 
a final state neutral kaon is identified in a Kg or a Kg state. This involves an interference 
between CF and DCS decays. A method for measuring this interference has been proposed 
in Ref. [36]. Defining a rate asymmetry between decays involving Kg and Kg, 


R(D°, M u ) = 


0 _ T(D° -> K° S M°) - T(D° -> K° l M°) 


T(D° -> KgM°) + T(D° ->• K° l M°) ’ 

one obtains, to leading order in the ratio of DCS and CF amplitudes, 

2Re [A(K 0 M°)A*(K 0 M°)\ 


(M° = 7r, r], Tj') 


R(D U , M°) = 


| A (/ 1 ° M °)| 2 


(79) 


(80) 


Eqs. (1761) predict equal asymmetries for M° = it and M° = r) 8 up to second order U-spin 
breaking. The two asymmetries are given by 


R(D°, V s) = R(D°, 7T°) = 2 tan 2 6 C 1 + 2Re + O(e 2 0 , e' 0 2 ) 


Comparing Eqs. (1781) and (ISTll we find 


|.4(iW)| _l r 

\A(K°M°)\ 2 R(D ,M )1 1 + °< € o. £ o) 


(M° = A m) ■ 


(81) 


(82) 


That is, although the branching fractions for DCS decays D° —> K°tt° and D° —>■ K 0r r]s 
cannot be measured directly, they may be obtained up to second order U-spin breaking 
corrections from corresponding CF branching fractions and Kg — Kg asymmetries: 


B(D° -> K°M°) = -[R(D°,M°)} 2 B(D° -> K°M °) [l + O(e 2 0 ,e' 2 ) 


(M 0 = tt,vs) • (83) 


Table I1VI summarizes current relevant information on branching ratios and amplitudes 
for D° decays into pairs of neutral pseudoscalars. We do not include the V6 mode (and 
decays involving the rf) , as D° —> rj 8 rj 8 would include D° —> r/r/ which has zero phase space. 
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Table IV: Branching fractions and amplitudes for D° decays to pairs of neutral pseudoscalar 
mesons 


Decay mode 

Branching fraction 

p (GeV/c) 

\A\ = \jB/p (GeV/c ) -1 / 2 

D° ->• A> u 

(1.19 ±0.04) x KT 2 

0.860 

— 

D° K° l tt° 

(1.00 ±0.07) x 10 ~ 2 

0.860 

— 

D° ->• /?°7T 0 

(2.29 ±0.07) x 10 “ 2 a 

0.860 

0.163 ±0.002 

D° -)■ K°r] 

(0.890 ±0.079) x 10 ~ 2 

0.772 

0.107 ±0.005 

D° 7r°7r° 

(8.20 ±0.35) x KT 4 

0.923 

0.0298 ± 0.0006 

D° —> rjir 0 

( 6.8 ±0.7) x 10 ~ 4 

0.846 

0.0284 ±0.0014 


“ Branching ratio calculated as twice the average of B(D° —> K^n 0 ) and B(D° —> K® 7 r°). 


An estimate of U-spin breaking is given by a ratio of SCS and CF decay amplitudes which 
equals one in the symmetry limit [see (R77|) (155]) ]. 


|A(7T 0 7r 0 )| 

|A(A'° 7 r°)| tan 6 c 


- 0.21 ± 0.02 , 


Another quantity measuring U-spin breaking is 


(84) 


V3\A(K°ri)\ 

\A(K°tt°)\ 


0.14 ±0.05 . 


(85) 


In order to determine Reeo from (1771) one would have to know also \A(K°r]')\ and the relative 
strong phase between this amplitude and A(K°r]), using a favored value for the 77 - 77 ' mixing 
angle [37] . The other U-spin breaking parameter, Re 6 q, is obtained from (| 8 TT) using a K G S —K° L 
asymmetry measurement by the CLEO collaboration [38], R(D°,n 0 ) = 0.108 ± 0.035, 


Regg = 


R(D°, 7T°) — 2 tan 2 6 c 
4 tan 2 6 c 


= 0.00 ±0.16 


( 86 ) 


Arguments favoring small U-spin breaking in the asymmetries R(D°, M°) have been pre¬ 
sented in Ref. [39] adopting a diagrammatic flavor SU(3) approach. 

The numerical values of Re e' 0 in (1861) and of the above other two measured U-spin breaking 
quantities imply that second order U-spin breaking terms in Eqs. (1771) - (l83l) are at most of 
order several percent. Neglecting these contributions involves an approximation that is 
quantitatively similar to the one used to obtain (|80l) . where terms which are second order in 
the ratio of DCS and CF amplitudes have been neglected. 


XII Conclusion 

We described a new approach to hadronic D° decay amplitudes applying a perturbative 
expansion in U-spin breaking parameters and treating isospin breaking carefully. We have 
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identified a class of two-body and quasi two-body decays involving charged pseudoscalars 
(P) and vector mesons (V), for which in each case adequate hierarchies have been shown to 
occur between U-spin breaking terms of increasing order. 

Nonlinear amplitude relations were predicted in each one of three cases, D° —> P + P~, 
D° —> V + P~ and D° —> P + V~ , which hold up to fourth order U-spin breaking and isospin 
breaking terms suppressed also by second order U-spin breaking. The three predicted rela¬ 
tions have been shown to hold experimentally at a very high precision varying between 1CU 3 
and 1CU 4 , in agreement with our estimates of high order terms. More precise tests require a 
first robust measurement of B(D° —>■ K + p~) and improving branching ratio measurements 
for D° decays to K* + n~, ir + K*~, n + p~ and K + K*~. So far no unexpected flavor symmetry 
breaking down to this very low level has been found, which would indicate physics beyond 
the standard model. This provides useful constraints on new | AC| = 1 operators, potentially 
originating in new physics at energies much above a TeV [TO] , 

Finally, we also studied decays to two neutral pseudoscalar mesons, deriving much less 
precise amplitude relations and relations for rate asymmetries between decays involving Kg 
and Kg, which hold up to second order U-spin breaking terms at a level of several percent. 

I wish to thank Yuval Grossman and Dean Robinson for useful communications. 
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